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Abstract 

Sudoku is a popular combinatorial puzzle. A new method of solving 
Sudoku is presented, which involves formulating a puzzle as a special type 
of transportation problem. This model allows one to solve puzzles with 
more than one solution, keeping the constraints of the problem fixed, and 
simply changing a cost matrix between solutions. 

1 Introduction 

Sudoku is a popular combinatorial puzzle consisting of a 9 x 9, partially- filled 
grid of digits 1 to 9. [1] The remainder of the grid must be filled in according 
to the rules: 

I. every row contains only one of each digit 1-9; 
II. every column contains only one of each digit; and 
III. every one of the nine 3x3 sub-grids contains only one of each digit. 

Sudoku has been modeled and solved using a wide variety of methods, for 
example, as an exact cover problem solved using the famous Dancing Links 
algorithm. [2] Examples from the field of operations research include the binary 
integer linear programQwhich involves 729 (81 x 9) binary variables representing 
the state of each of the 81 places in the square. [5] This is equivalent to a 
constraint satisfaction problem, where any feasible solution is a solution. [3] 

The latter methods are equivalent to the general transportation problem. 
However, for problems with more than one solution, these methods require 
iteratively reformulating the problem by inserting additional constraints until 
no new solution is found. The method proposed here differs from these in that 
the transportation problem provides an intuitive and compact representation 
for solving problems with many solutions. This method distinguishes between 
different solutions to a given problem by manipulating the cost matrix, without 
otherwise modifying the constraints of the problem. 
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Not to be confused with "binary Sudoku" 



http : //xked . com/74/ 



2 Sudoku as a transportation problem 



2.1 A general transportation model 

Transportation problems are a special type of linear programming problem. 
Ch. 14, p. 299] A transportation problem generally consists of: 

• a set of to supply points, s = {s^; i = 1, . . . , m}; 

• a set of n demand points, d — {dj] j = 1, . . . ,n}; and 

• a variable cost matrix, c = {cjj-; « = 1, . . . , m, j = 1, . . . , n}, representing 
the cost of transporting one unit from each supply point to every demand 
point. 

The objective is to minimize the total cost of shipments while meeting de- 
mand. Let Xij be the number of units shipped from a supply point i to a 
demand point j. Then a transportation problem has the general formulation: 
§7-1] 

m n 

minimize z = CijXij (1) 

i=i j=i 



subject to 2_j Xi 'i — Si Vi = 1, 
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^ Xi,j >dj Vj = l,...,7i 

i=l 

Xij > Vi=l,...,m;j = l,...,n 

If the total supply of goods is equal to the total demand, that is Sj = 
J2j dj, then the problem is said to be balanced. 

Solving a Sudoku puzzle involves placing (or transporting) up to eighty one 
digits into an equal number of cells. Therefore, in modeling a Sudoku puzzle 
as a transportation problem, there are eighty one demand points to consider, 
each with a demand of one unit. The eighty one digits available to be placed 
in the grid consist of (nine of each of) the digits one through nine. Thus, there 
are nine units in each of nine distinct supply points to be considered. That is, 
m = 9 and n = 81 in Equations [T] Let s be the set of supply points, and let 
d be the set of demand points, such that s$ = 9 for all i = 1, . . . , 9 and dj = 1 
for all j = 1, ... ,81. An illustration of this transportation network model of 
Sudoku is shown in Figure |2.1| 



2.2 Transportation cost matrix 

In modeling different Sudoku puzzles, the supply and demand sets s and d 
remain constant. A variable cost, Cij, is incurred when transporting one unit 



2 



di = 1 



dj 



81 

d 81 = 1 



Figure 1: A transportation network model of a Sudoku puzzle. There are nine 
supply points, with a supply of nine units each; and eighty one demand points, 
with a demand of one unit each. Each one of the nine supply points, i, can ship 
to every one of the eighty one demand points, j, one unit at a cost of Cij. 



from a supply point i to a demand point j. Thus, it is necessary to formulate 
a cost matrix c which assigns a cost of transporting one unit from supply point 
to demand point. 

Let a puzzle matrix, II, be a 9 x 9 matrix representing the contents of 
the initial Sudoku puzzle grid. An element of the matrix II is non-zero if its 
corresponding cell in the Sudoku grid contains a digit, in which case its value is 
equal to that of the digit in the cell, and zero otherwise. Let 7Tjj represent the 
elements of the puzzle matrix II, where i = 1, . . . , 9 and j — 1, ... ,9. 



n = 



TTi l • • • 7T1 



7I"9,1 ' ' ' To 9 



Each element in this matrix has an integer value A = 0,1,..., 9. Given 
an initial puzzle matrix II, we seek Cij for all i = 1, . . . , 9 and j = 1, . . . , 81. 
Consider one element in a given puzzle matrix, ir a ^ = A. If A is non-zero, 
then according to rules I and II, the following elements of II cannot then also 
have a value of A: all elements (a, j), j = 1, . . . , 9 ^ 6; and all elements («,&), 
i = 1, . . . , 9 ^ a. Furthermore, consider that for any coordinate pair (a, b) where 
a = 1, . . . , 9 and b — 1, . . . , 9, there exist p, q, r, and s, such that a — 3p + q, 
where p = L a /3J, q = a mod 3; and b = 3r + s, where r = [b/3\, s = b 
mod 3. Thus, according to rule III, A cannot also be the value of any element 
(3p + to, 3r + n) for all m = 1, 2, 3 and n — 1,2, 3, where p = [a/3\ , r = [b/3\ , 
and where 3p + m ^ a and 3r + n ^ b. 

Thus, given some non-zero element in II, 7r a ^ = A ^ 0, let R (a, b) be defined 
as the set of all other elements in II whose value cannot also be equal to A, as 
determined by rules I to III: 
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R(a,b) = {(i,b)A*J), (3 UJ 



m, 3 



+ « \{(M)} (2) 



where i — 1,...,9, j = 1,...,9, m = 1,2,3, and n = 1,2,3. This set 
represents the set of all cells in the same row and column, and all cells within 
the same 3x3 sub-matrix, as Tr a ^, excluding itself. Let this set be called the 
restricted set. 

An illustration of the restricted set for a single digit in a Sudoku grid is 
shown below; in this example, ^4.4 = 3 and the corresponding R (4, 4) is shown 
on the right. 

There exist as many restricted 
sets as there are non-zero elements 
in II. The intersection of all the re- 
stricted sets for any given II can be 
expressed in the form of a matrix, 
which will be used to formulate the 
cost matrix. 

Consider all non-zero elements 
of the puzzle matrix II, 7T 0l b = A 7^ 
0. For each of these, let cj)^ a ^ be a 
matrix of dimensions 9x9, defined as: 



-M 
= { M 



if {hi) = (a,b) ; 

if (i, j) £ R (a, b) ; and 

otherwise; 



(3) 



where i = 1, . . . , 9, j = 1, . . . , 9, and M > 0. 

Consider the nine possible values, A = 1, ... ,9, of each non-zero element of 
II. For each, let cj)^ x ' be a matrix of dimensions 9x9, defined as: 



(a,fc):7r a b— A 



Finally, let represent the row vectors of 4>^ x \ where i 
<I>^ be the vector defined as: 

*W = [n ••• r 9 ]. 
The variable cost matrix sought, c, is given by: 

$(i) 



$0) 



(4) 

, 9, and let 
(5) 

(6) 



Consider again Equation [3] A cost of M means that it is highly undesirable 
to ship a unit to a particular demand point. Conversely, a cost of — M means 
that the shipment is highly desirable. A cost of zero means that the shipment is 
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possible but without preference over other possible destinations that also incur 
a zero cost. To understand the meaning of M intuitively, consider the limit 
M — » oo. If Cjj = M, this means supply point i can not ship to demand point j 
"at any cost." Conversely, if Cjj = —M, this means supply point i must ship to 
demand point j "at all cost." In practice, the value of M need only be strictly 
positive. Algorithms that solve transportation problems are guaranteed to find 
an optimal solution if one exists. An optimal solution will not assign a shipment 
from a supply to a demand point when the cost of such a shipment is M; such 
shipments will all be avoided. In fact, any optimal solution will only contain 
shipments between points that incur zero cost. For simplicity, let M = 1. 

There can be more than one optimal solution to the general transporta- 
tion problem posed by the supply vector s, the demand vector d, and the cost 
matrix c. These solutions all have an objective function value of z = —NM, 
where N is the number of digits initially given in the Sudoku puzzle. However, 
these solutions do not all form Latin squares, much less valid solutions to the 
Sudoku puzzle. If the Sudoku puzzle has a unique solution, then exactly one 
of the optimal solutions to the general transportation problem will also be the 
solution to the puzzle. Therefore, additional constraints to those in the general 
transportation problem are needed. 

2.3 A special transportation problem 

We formulate rules I to III as additional constraints, to obtain the special trans- 
portation problem: 



minimize z — ^ ^ c i,j x i,j (?) 
i=i j=i 

n 

subject to ^ Xij — Si i = 1, . . . , m 



^ ^ dj j 1, . . . , Tt 

i=l 
n 

^ Xij < 9 i = 1, . . . ,m 

i=i 

m 

^2xij<l j = l,...,n 

i=l 

Xij G {0, 1} i = l,...,m; j = 1,... ,n 

It is now possible to solve the Sudoku puzzle by solving the transportation 
problem posed by the supply vector s, the demand vector d, and the cost matrix 
c as defined above, and these additional constraints, using any one of well-known 
algorithms, such as the simplex algorithm. This is an efficient way of solving 
Sudoku puzzles. 
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3 Conclusions 



It has been shown that a Sudoku puzzle can be modeled as a special type of 
transportation problem. The model allows to efficiently solve for all solutions 
to any given problem, by keeping constraints fixed and manipulating the cost 
matrix. Algorithms for solving transportation problems are well established and 
have the attractive property of generally completing in polynomial time. 

Sudoku is also a special case of a greater class of problems - that of com- 
pleting partially-filled Latin squares - which has applications in cryptography. 
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